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ABSTRACT 

We argue that the bright flare of the binary pulsar PSR B1259— 63/LS2883 
^ ■ detected by the Fermi Large Area Telescope (LAT), is due to the inverse Compton 

(IC) scattering of the unshocked electron-positron pulsar wind with a Lorentz 
factor To ~ 10 4 . The combination of two effects both linked to the circumstellar 
disk (CD), is a key element in the proposed model. The first effect is related to 
the impact of the surrounding medium on the termination of the pulsar wind. 
Inside the disk, the "early" termination of the wind results in suppression of 
its gamma-ray luminosity. When the pulsar escapes the disk, the conditions 
for termination of the wind undergo significant changes. This would lead to a 
dramatic increase of the pulsar wind zone, and thus to the proportional increase 
of the gamma-ray flux. On the other hand, if the parts of the CD disturbed 
by the pulsar can supply infrared photons of density high enough for efficient 
Comptonization of the wind, almost the entire kinetic energy of the pulsar wind 
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would be converted to radiation, thus the gamma-ray luminosity of the wind 
could approach to the level of the pulsar's spin-down luminosity as reported by 
the Fermi collaboration. 



Subject headings: binaries: close - - gamma rays: stars - 
nisms: non-thermal — pulsars: individual (PSR B1259— 63) 
(LS2883) 



radiation mecha- 
- stars: individual 



Introduction 



The pulsar magnetospheres are effective gamma-ray emitters (jAbdo et al.1 I201CH ) . Pul- 
sars can produce potentially detectable gamma-ray emission also due to the bulk Comp- 
tonization of their cold ultrarelativistic winds. While in the case of isolated pulsars this 
radiation c omponent is generally weak, unless the wind is accel erated (relatively) close to 
the pulsar (JBogovalov k Aharonianl l2000l : lAharonian et al.l 120121 ) , in binary systems the ra- 
diation is significantly enhanced due to the presence of dense target photon field supplied b y 
the optical companion ( ball k Kiri3l200ot ball k Pod JiooJ JKhangulyan et al.ll2007l . l201 if ). 
Gamma-rays are produced also after termination of the wind due to the IC scattering of rel- 
ativists electrons accelerated by the termination shock. In isolated pulsars the termination 
of the wind leads to the formation of the so-called Pulsar Wind Nebulae with a distinct 
extended synchrotron and IC emission which can be easily separated from the radiation of 
the unshocked wind. In binary pulsars, the separation of the radiation components produced 
before and after the termination of the pulsar wind is a more difficult task; it requires careful 
analysis of the spectral and temporal features of two radiation components. 

The most promising object for exploration of processes of formation, acceleration and 
termination of pulsar winds in binary systems is PSR B1259— 63/LS2883. It contains a 
47.7 ms pulsar orbiting a luminous star in a ver y eccentric orbit (eccentricity e = 0.87, period 
-Porb = 1237 d, semi-major axis a 2 = 7.2 AU, see lJohnston et al.lll992uNegueruela et al.ll2011 



and references therein). Th e variable TeV gamma - ray em ission detected from this object by 
the H.E.S.S. collaboration (Aharonian et al.ll2005l 120091 ) was in general agreement with the 
predictions fjKirk et al.lll999i ). but the TeV gamma-ray lightcurve appeared to be different 
from expectatio ns. Several possible scen arios suggested for interpretation of the reported 
TeV lightcurve (jKhangulyan et al.l 120071 ). predict essentially different gamma-ray fluxes at 
MeV/GeV energies, therefore the measurements in this energy band should significantly 
contribute to the understanding of the hydrodynamics and the particle acceleration processes 
in binary pulsars. 
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The recent observations of PSR B1259— 6 3/LS2883 with Fermi LAT revealed a complex 
behavior of the system in GeV gamma-rays (jAbdo et al.ll201ll ; iTam et al.ll201ll ). Close to 
the periastron passage, the source has been detected with a modest energy flux of 6 x 
10~ n ergcm~ 2 s _1 . However, +30 days after periastron passage (with largest recorded flux 
level at +35 days), a spectacular flare has been recorded. It lasted approximatel y two weeks 



with an enhanced flux above 100 MeV at the level of 3 x 10 10 ergcm 2 s l ( lAbdo et al. 



201 ll ; ITam et al.l 120111 ) . The flare is characterized by a sharp increase and a rather smooth 
decay over approximately 2 weeks. Importantly, the strong i ncrease of gamma ray flux was 
not accompanied by a noticeable change of the X-ray flux ( lAbdo et al.l l201lh . While the 
weak GeV flux reported during the periastron passage can be explained by IC scattering of 
electrons accelerated at the termi nation shock and/or by the bulk Comptonization of the 
unshocked electron-positron wind (JKhangulyan et al.l 1201 if ), the post-periastron flare was a 
real surprise. In general, this flare represents a unique case in astrophysics when the available 
energy of an object is fully converted to nonthermal high energy radiation. On the other 
hand, given the current belief that the rotational energy of the pulsar is converted to kinetic 
energy of cold ultrarelativistic electron-positron wind, the binary pulsars can in principle 
operate as perfect '100 %' efficient gamma-ray emitters if the density of the surrounding 
target photon field would be high enough for realization of Comptonization of the pulsar 
wind in the saturation regime. Yet, any explanation of this flare should address several 
important issues, namely (i) the extraordinary high luminosity of gamma radiation at the 
level close to the pulsar's spin-down luminosity (SDL), E = 8 x 10 35 erg/s; (ii) lack of 
enhanced non-thermal X-ray flux; (iii) orbital phase of the flare. 



The scenario 



The very fact of absence of simultaneous activity in the X-ray energy band suggests 
that the gamma-ray emission is produced by the pulsar wind. Indeed, since in the cold wind 
the electrons move together with magnetic field, they lose their energy only through the IC 
scattering, thus IC gamma-rays are emitted without accompanying synchrotron radiation. 
Moreover, for a narrow, e.g. Planckian type distribution of ambient photons with temper- 
ature T r , the upscattered gamma-ray emission will be concentrated within a narrow energy 
interval with a characteristic energy E 1 ps 3kT r Tl « 300(A;T r /l eV)TQ 4 MeV in the Thom- 
son regime (A;T r r /m 2 c 4 <c 1) or E 1 s=s m e c 2 T « 5r 0i 4 GeV in the Klein-Nishina regime 
(A;T r r /m 2 c 4 > 1). Here T = 10 4 r 04 is wind bulk Lorentz factor. 

The strength of the gamma-ray signal produced by the wind depends on three pa- 
rameters: (i) the bulk Lorentz factor of the wind, (ii) the density of the target photon 



-4- 



field, (iii) the pulsar wind zone (PWZ) length towards the observer ( s ee e.g. iBall &: Dodd 



2001 



Khangulvan et a 



2007t ISierpowska-Bartosik fc Bednarekl l2008t ICerutti et al 



2008 



Khangulvan et al.ll201ll ). At least one of these parameters should experience sudden changes 
in order to prov i de a sharp increase of the gamma-ray signal as detected by Fermi LAT. 
Bogovalov et al.l (120081 ) have shown that the interaction of the pulsar wind can proceed in 
two different regimes depending on the ratio 77 of ram pressures of the interacting pulsar 
and stellar winds. For rj < r] CT , where r] cr ~ 10~ 2 , the pulsar wind termination shock has a 
closed structure, and for r] > r] rr the termination shock is expected to be unclosed (on a scale 
of the binary system, see iBosch-Ramon fc Barkovl l201ll ; iBosch-Ramon et al.l 120121 ) . This 
hydrodynamical instability should lead to a fast (on timescales of 10 3 — 10 4 s) transformation 
of the termination shock. Importantly, in the case of interaction of the pulsar with the CD 
r] d rsj 10~ 3 <C r] CT , w hile at the interaction o f the pulsar wind with the stellar wind (SW) 
i| ff ~5x 10~ 2 ^> r] cr (jKhangulyan et al.ll201ll ). Thus, the pulsar's entrance to and exit from 
the CD should lead to fast transformations of the termination shock. In Figure HJ we show 
PWZ lengths calculated for different orbital inclinations and the rj parameter (for details see 
Khangulvan et al.l 1201 ll ). 



The increase of the PWZ length by a factor of ~ 10 after its escape of the disk naturally 
explains the detected dramatic increase of the gamma-ray flux during the flare (see Fig. [2]). 
What concerns the absolute gamma-ray flux with corresponding luminosity comparable to 
the pulsar's SDL, it requires that the Comptonization of the wind electrons proceeds in the 
optically thick regime, i.e. the IC cooling length lie should not exceed the PWZ length of 
about lie ~ 10 13 cm outside the disk. In the Thomson regime, 



he ^ 10 14 r ^(w ph /ergcm 3 ) W. 



T 



where w p h is target photons energy density. Thus, the required energy density of the target 
photon field is 

w ph ~ 10T 0i ll~ c \ 13 erg cm" 3 , (2) 

where hc,i3 = ^Tc/(10 13 cm). Althou gh the luminosity of the optical companion LS 2883 is 
rather high (iNegueruela et al.l 120111). it is s till not sufficient for realization of the optically 
thick IC scenario (jKhangulyan et al.l 120111). The CD is another source of target photons 
for Compton scattering of electrons (Ivan Soelen fc Meintjesl l201ll ) , which contribution can 
be significant (especially in the close vicinity of the disk). Equations flTJ and (j2|) give an 
estimate of the required luminosity of the target photon field 



L ph >10 38 r^/ IC ,i 3 ergs- 



(3) 



This is a quite tough req uirement, as large as 40 % of the luminosity of the optical companion 
(INegueruela et al.l l201ll ) . One may speculate that such a high photon field is a result of 
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thermal radiation of the heated parts of the CD triggered e.g. by the interaction of the 
pulsar with the disk. Whether such an effective heating of the disk is possible is an open 
issue. While the answer to this question requires detailed hydro dynamical studies, below we 
simply postulate the existence of such a high density radiation field. 

The requirement on the luminosity of the target photon field given by Eq. fl3]), can be 
relaxed assuming larger wind bulk Lorentz factors and/or smaller values of /k> The Lorentz 
factor cannot much exceed 10 4 , unless the temperature of radiation is significantly less than 
10 4 K (to explain the reported energy spectrum of gamma-rays). However, for T r <C 10 4 K, 
the required energy density of the photon field would exceed the Stefan-Boltzmann limit for 
black-body radiation. Regarding the size of the region filled with photon field, it cannot 
be smaller than > 10 13 cm, given that during the flare (from +30d to +60d from periastron 
passage) the star separation distance changes from 3.7 x 10 13 cm to 6 x 10 13 cm. 

In the framework of the proposed model, one should expect a similar flare when the 
pulsar exits the disk during the pre-periastron epoch. However, due to the differences in 
the termination shock orientation, the intensity of this flare should be different (see Fig 12]). 
Assuming that the pulsar exits the disk at epochs symmetric in terms of true anomaly, the 
pre-periastron exit should occur 10 days before the periastron passage. For rj = 5 x 10 -2 the 
ratio of PWZ lengths corresponding to +30 d and — 10 d epochs exceeds 4 (see Figs. [Hand 
[2]), therefore the expected flux level is consistent with the Fermi LAT upper limits obtained 
for that period. 

Below we present and discuss the results of numerical calculations performed within the 
suggested scenari o. The formalism appli e d to PSR B1259— 63/LS2883 is described in our 
previous papers (jKhangulyan et al.l 120071 . 1201 1[ ). The calculations rely on h ydrodynamical 
simulations of interaction of the pulsar and CD (jBogovalov et al.ll2008l . l2012h . In this paper 
we include the effects related to the presence of the CD which have not been taken into 
account in previous studies. 



3. Radiation of the unshocked pulsar wind 

In Figure H] we show the spectral energy distribution (SEP) of IC radiation of the 
pulsar wind averaged over the period of Fermi LAT observations (jAbdo et al.ll201ll ). The 
calculations are performed for different values of the initial Lorentz factor r and the rj - 
parameter. These value correspond to the following regimes: rj — 1 - up per limit; r? = 0.0 5 



- int eraction with the SW; rj = 10 3 - interaction with the CD (see iKhangulyan et al. 



20 111 , for details). The comparison of calculations with the Fermi data shows that the IC 
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radiation of the unshocked wind can achieve the level of reported fluxes. However, these 
fluxes are also rather close to theoretical predictions for radiation of the shocked wind (see 



Figs. 7 and 13 in iKhangulyan et al.l 120071 ). Unfortunately the marginal gamma-ray signal 
obtained during the periastron passage does not allow a quantitative analysis which would 
give preference to the radiation associated with the wind before or after its termination. 
On the other hand, the detected gamma-ray fluxes can be treated as upper limits for the 
unshocked wind component. This allow us to constrain the wind's Lorentz factor T ; as it 
follows from Figure HJ the interval of r ~ (0.1 — 3) x 10 5 is excluded by the gamma-ray 
data. 

The gamma-ray signal from the pulsar wind associat ed with the target phot ons of the 
companion star has a smooth orbital phase dependence (IKhangulyan et al.l 1201 ll ). In con- 
trast, the CD can introduce quite sharp tempora l features. Presently , the properties of CDs 
in binary pulsars are a subject of debate (see e.g. lOkazaki et al.ll201l[ ). Here we assume that 
a CD does exist in PSR B1259-63/LS2883, and that it extends to 3 x 10 13 cm from the 
star, as it follows from the eclipse of the pulsed radio emiss ion observed in 2010 between 
—16 and +15 days around the periastron (jAbdo et al.ll201ll ). The impact of the disk has 
been taken into account by introducing two effects: (i) the change of the ^-parameter; and 
(ii) the presence of an additional target photon field from the disk itself. 

To illustrate the importance of these effects, we performed calculations corresponding 
to the three important cases: (a) suppression of the PWZ length by the disk, i.e. rj = 10~ 3 , 
t = +35 days after periastron passage; (b)flare epoch, i.e. 77 = 0.05, t = +35; and (c) 
pre-periastron flare, i.e. r] = 0.05, t = —10 days before periastron passage. The additional 
target photon field was assumed to be gray body with temperature T r = 10 4 K and energy 
density u> p h = 2.8ergcm~ 3 . The results of calculations are shown in Fig. [2j The sudden 
increase of the gamma-ray flux by an order of magnitude during the flare is due to the 
increase of PWZ length just after the pulsar leaves the disk. During the first 1-2 weeks of 
this epoch the pulsar is located close to the disk, and therefore the wind zone is intensively 
illuminated by the photons of the disk. Apparently, the dominance of the infrared component 
of the disk over the optical radiation of the star, and the upper limit on the Lorentz factor 
of the wind derived from Fermi LAT observations during the periastron passage, explain 
uite naturally th e significantly sof ter spectrum of gamma rays observed during the flare 
Abdo et al.ll201ll ; iTam et al.l 120 111 ). From Fig. [2] one can see that with the assumed set 
of three model parameters rj, T r , and To, one can explain both the spectral shape and the 
absolute flux of gamma-rays during the flare. With a slightly different choice of parameters 
than the ones used in Fig. [5] we can provide, in principle, a perfect fit to the observational 
results. However, given the uncertainties related the complex character of the interaction 
of the wind with the disk, the "perfect fit" would be redundant. Instead, we would like to 
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emphasize the most important consequences of this interpretation. First, it implies that the 
Lorentz factor of the wind cannot be significantly smaller than 10 4 given that the energy 
of the IC photons scales as E 1 ex T v Tq. On the other hand, the Fermi LAT observations 
during the periastron passage constrain the wind Lorentz factor to r < 10 4 . Thus, if our 
interpretation of the flare is correct, the wind Lorentz factor should be quite close to 10 4 . 



4. Discussion 

In this letter we propose a model for the recent Fermi observations of gamma-rays from 
PSR B1259— 63/LS 2883. We assume that gamma-rays are produced at interaction of the 
cold pulsar wind with external photon fields. While the Comptonization of the wind by 
radiation of the companion star can explain the modest gamma-ray flux detected during the 
periastron passage, the order of magnitude higher flux of the post-periastron flare requires an 
additional seed photon component. We propose that the radiation of the heated CD can play 
this role. The Comptonization of the pulsar wind by the CD emission proceeds in a quite 
specific way. Inside the disk it is suppressed because of the high ram pressure in the disk. 
The density of the stellar optical photons generally dominates over the density of IR photons. 
But in the immediate vicinity of the disk, where the IR density can be still very high, the 
optical depth for Compton scattering might be significantly larger than 1. This would lead 
to a very effective, close to 100% efficiency of transformation of the rotational energy of the 
pulsar to high energy gamma rays via Comptonization of the cold ultrarelativistic wind. In 
the suggested scenario, the departure of the pulsar wind from the disk will lead to a gradual 
decrease of the photon density from the disk and correspondingly to the reduction of the IC 
gamma-ray flux on time-scales of days. 

During the exit of the pulsar from the disk at the epoch before the periastron passage, 
the termination shock is expected to expand towards the direction opposite to the observer, 
thus we should not expect strong gamma-ray signal. This agrees with the Fermi LAT 
observations. The main issue for realization of the proposed scenario is whether the pulsar- 
disk interaction can provide the required photon field. This is a key question which should 
be answered by detailed studies of interaction of the pulsar wind with the CD. This study 
should also address another important issue, which is related to the 15 days delay between 
the appearance of the pulsed radio emission and the epoch of the GeV flare. 

If the proposed interpretation is correct, this would be the second case of direct mea- 
surement of a pulsar wind's Lorentz factor. Recently, the Lorentz factor of order of 10 6 has 
been extracted fr om the observation of very high energy pulsed gamma-ray emission from 



the Crab Pulsar (jAharonian et al.ll2012l ). Interestingly, the two orders of magnitude smaller 



Lorentz factor of the pulsar wind in PSR B1259— 63/LS 2883 compared to the Crab pulsar, is 
quite close the ratio of spin-down luminosities of two pulsars. If this is not just a coincidence, 
the linear dependence T oc L$n would imply a universal production rate of electron-positron 
pairs of order of 10 38 e ± /sec, independent of the pulsar spin-down luminosity. 

Finally, we should note that in principle the gamma-ray flare can be explained also by 
the post-shock flow. Since the flux detected by Fermi LAT was close to pulsar's SDL, this 
interpretation requires significant Doppler boosting o f radiation. Relativist ic flows can be 
indeed formed at interactions of the pulsar and SWs (JBogovalov et al.l 120081) . and therefore 



the broad-band ra diation can be affected by the Doppler boosting (IKhangulyan et al.l 12008 



Dubus et al.ll2010l ). The Doppler boosting should amplify also the X- and TeV gamma ray 



fluxes. However, the lack of any noticeable activity during the flare at other wavelengths 



(lAbdo et al.l 120 111 ) makes this interpretation less likely. Also, one should note that the flux 
level of the Doppler boosted radiation can easily exceed the limit set by the pulsar's SDL, 
but flares with apparent super SDL so far have not been detected. On the other hand, the 
detection of such events in future observations would rule out the origin of gamma-radiation 
related to the unshocked pulsar wind. 
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Fig. 1. — Upper panel: PWZ lengths towards the observer as a function of days to periastron. 
The calculations are performed for two values of the ^-parameter: r\ = 10~ 3 (CD) and 
7] = 0.05 (SW) and three different orbital inclinations: % = 25°, i = 45°, and i = 75°. 
Bottom panel: ratio of the PWZ lengths for the interaction with SW and CD for —10 d 
(dot-dashed line) and +30 d (dash-dot-dotted line) to periastron passage; and expected 
ratio of the post- to pre-periastron flares (dotted line). The inset plot shows dependence of 
the true anomaly (in degrees) on the epoch (in days) to periastron. 



JO 



C\J 



I 10" 10 

D) 


x" 



10 



11 



10 



- 12- 



10" 3 - 



12 



r|=10" 3 ;dt=+35d ■ ■ 
ti=0.05; dt=+35d — 


i | i i 


Abdoetal. (2011) • 


"r(=0.05;dt=-10d 




5 


; r 0= 10 4 


T =4.6x1 o 4 


r =2.2xi o 5 *j£--n - 

.^^ v X^ * 1 

/^ ' > X ' •! 

X * \v * * I 




X «»""5 


&/ / .* / .* \ 


s h * H 


Y ■'? 


v / -" / V / .* 1 

♦ \ \ / ♦ X 1* * • 


7 / '"} 


r\ S ♦ / 

Yt * .* x 


*\ X' ♦* / U * ♦* 

A ' • •• X 1. ' ♦ 

V y / w / 




/ * .* / 


• ' f > v X 


/i .i x K .♦ 




/ * ♦ / 

/// 
// / 


♦ V ♦ 

/ ♦ • \f 

* *r 

♦ «n 

♦* XM 


•\ X 

1 1 ^. ' . 


♦1 X 

• ■ x B 
• lx ' 

♦* X ■ 

XI ' 


* ♦i | 
^ ♦ r. 


- 


1 > ' 
1 

1 4 
1 ' 


/ ♦ 

* 
♦ i 

♦ X 
* X 


/ I / ♦ 


II ! - 




/ ♦ 


■ 




. / * ♦ i 




/• ♦ 






,/ » 






* / ' .* x 




/ »* 




1 


/ 


■ 


_ 


/ * ' / 




* • i 




1 


1 . 


■ 




r / .J/J 


4 
1 

. /I 


•* •> 

♦ «X 




, A- 


« li 


\ ■ " 



0.1 



10 



100 



1000 



E yJ GeV 



Fig. 2.— The spectral energy distribution of the IC radiation of the unshocked pulsar wind 
at the epoch of 35 days after periastron passage [r\ = lCT 3 -dotted lines; and r\ = 0.05-solid 
lines) and at the epoch of 10 days before periastron passage (rj = 0.05 - dashed lines). 
The calculations are performed for two target photon fields: (i) radiation of the companion 
optical star with a radius i?* = 6.2 x 10 11 cm and the surface temperature T, = 3x 10 4 K, 
and (ii) radiation of the CD, which was assumed to be isotropic gray body with temperature 
T r = 10 4 K and energy density w p h = 2.8 erg cm -3 . Several initial pulsar wind bulk Lorentz 
factors have been assumed: T = 10 4 , 4.6 x 10 4 , 2.2 x 10 5 and 10 6 . It is assumed that 
the kinetic energy luminosity of the cold ultrarelativistic pulsar wind is equal to the SDL 
of the pulsar. The calculations performed for the model parameters rj = 5 x 10~ 2 and 
To = 10 4 match quite well the observational points ( pentagons) report ed by the Fermi LAT 
collaboration for the post-periastron flaring episode (JAbdo et al.ll201ll ). 
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Fig. 3.— Sketch of the scenario. 
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Fig. 4. — Spectral energy distributions of IC radiation of the unshocked pulsar wind averaged 
over the period from —35 to +18 days around the periastron passage. The calculations are 
performed for different values of the ^-parameter: r\ — 1 (solid lines), r\ = 0.05 (dotted 
lines, SW) and r\ = 10 -3 (dashed lines, CD), and the pulsar's wind bulk Lorentz factor: 
T = 10 4 , 3 x 10 4 , 10 5 and 3 x 10 5 . The densities of target photons at different positions of 
the wind are calculated assuming a spheric al star with a radius R* = 6.2 x 10 11 cm and surfac e 
temperature T* = 3 x 10 4 K (see for details iNegueruela et al.ll2011uKhangulyan et al.ll201ll ). 
Th e gamma-ray meas urements correspond to the period of [—35, 0] days (pentago ns) reported 



bv iTam et al.l (]201ll ). and the period of [—28,18] days (squares) reported by lAbdo et al. 
J201 il l. 



